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Abstract 

The mechanism of the Coulomb breakup reactions of the nuclei with neutron- 
halo structure is investigated in detail. A time-dependent Schrodinger equa- 
tion for the halo neutron is numerically solved by treating the Coulomb field 
of a target as an external field. The momentum distribution and the post- 
acceleration effect of the final fragments are discussed in a fully quantum 
mechanical way to clarify the limitation of the intuitive picture based on the 
classical mechanics. The theory is applied to the Coulomb breakup reaction 
of ^^Be -|- ^'^'^Pb. The breakup mechanism is found to be different between the 
channels of j'" = ^ and | , reflecting the underlying structure of ^^Be. The 
calculated result reproduces the energy spectrum of the breakup fragments 
reasonably well, but explains only about a half of the observed longitudinal 
momentum difference. 

PACS number(s): 25.60.-Fv, 25.70.De, 25.70.Mn 
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I. INTRODUCTION 



The neutron-halo structure has been observed systematically [|T| in light neutron-rich 
nuclei near the neutron drip-line. A large Coulomb breakup cross section has been observed 
for the neutron-halo nucleus in reactions on a heavy target nucleus , which indicates that a 
significant amount of El strength exists at low excitation energy region of the neutron-halo 
nucleus [Q. Since stable nuclei do not have such a strong El distribution at low excitation 
energy, this unusual feature of the neutron-halo nucleus has attracted much attention. The 
origin of the strong El distribution is still controversial: It may be a resonant character due 
to the vibration of the halo neutron against a core nucleus . Or it may be understood from 
an analogy of the strong El transitions known in ^Be, ^^Be, and ^^C where a single nucleon 
is weakly coupled to a core nucleus A recent argument suggests that the low-lying El 
strength of the light halo nuclei is not considered the vibrational state 

Recent experiments of the Coulomb breakup reactions, ^^Li -|- ^°^Pb and ^^Be +^°^Pb 
P], have observed a significant longitudinal momentum difference between the halo neu- 
tron(s) and the core nucleus. The momentum difference has been explained in terms of 
the Coulomb post-acceleration effect by assuming a direct breakup mechanism. In this 
mechanism the breakup is assumed to occur instantaneously at the closest approach point 
between the projectile and target nuclei. After the breakup of the projectile nucleus, the 
target Coulomb field accelerates only the core nucleus, and causes the momentum difference 
between the neutron and the core nucleus. If the breakup proceeds through a resonant 
state of the projectile nucleus, the core nucleus and the halo neutron moves together during 
its lifetime and the post-acceleration effect should become small. The observation of the 
longitudinal velocity difference is thus recognized as a direct evidence for the non-resonant 
character of the El strength at low excitation energy. 

The post-acceleration effect is not explained in the lowest order perturbation treatment 
of the Coulomb excitation. The evaluation of higher order effects is not easy because the 
final states involve continuum states. Several theoretical approaches have been proposed 
to understand the post-acceleration effect, including a classical treatment of the breakup 
reaction 0, a distorted- wave Born approach |jTO|, a simplified treatment of the higher order 
perturbations [O], and a coupled-channel approach with discretized continuum states [IT2 



Contrary to these approaches, some groups have investigated the time evolution of the 
projectile nucleus by solving a time-dependent Schrodinger equation on mesh points of space 
and time variables |T3| , p^ . 

In our previous paper |]14[ we studied the Coulomb breakup of ^^Be and found large 



transverse and small longitudinal momentum differences between the neutron and the ^°Be 
nucleus. The result was understood in the picture of free-particle breakup mechanism which 
is in contrast to the direct breakup mechanism. However, the reproduction of the experi- 
mental momentum difference remained to be an open problem. In this paper we extend the 
previous calculation to a more realistic case by including the spin-orbit interaction between 
the neutron and the core nucleus. A full three-dimensional dynamical calculation has also 
been done in [13b] assuming a simple internal Hamiltonian. No investigation has, however, 
so far been performed to clarify the roles of the Is potential and the level structure of excited 
states. We first analyze the mechanism of the Coulomb breakup quantum-mechanically for 
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various cases of the neutron-core Hamiltonian. The results are discussed in comparison with 
the intuitive arguments based on the classical mechanics. The usefulness and limitation of 
the classical arguments is made clear. We then analyze the Coulomb breakup reaction of 
^^Be + ^'^^Pb PI with a realistic choice of the potential between the halo neutron and ^°Be. 
We show that the breakup mechanism is sensitive to the structure of the excited states of the 
projectile and that the inclusion of the spin-orbit interaction is very important to reproduce 
the observed features quantitatively. 

In section 2 we formulate the quantum-mechanical treatment of the Coulomb breakup 
reaction using the time- dependent Schrodinger equation. The method of calculating the 
time evolution of the wave function is briefly explained. In section 3 the theory is applied 
to various cases of the halo nucleus to discuss the limitation of the classical arguments for 
the Coulomb breakup reaction and to reveal the characteristics of the quantum aspect. The 
Coulomb breakup of ^^Be on a ^°'^Pb target is analyzed in section 4. A brief summary is 
given in section 5. 



II. FORMULATION 



We consider the Coulomb breakup reaction of the nucleus with neutron-halo structure. 
The projectile of the halo nucleus is assumed to consist of a single neutron and a core 
nucleus. The core nucleus is treated as a structureless particle and binds the neutron weakly 
by an appropriate potential. We describe the reaction in the projectile rest frame where the 
center-of-mass of the projectile is put at the origin of the coordinate. The time development 
of the wave function, \E'(r, t), of the relative motion between the neutron and the core nucleus 
is described by the following time-dependent Schrodinger equation, 

zh^^^{r,t) = {if(r) + Kxt(r,t)}v|/(r,t), (1) 

where H(r) = — |^V^ + V{r) is the internal Hamiltonian describing the relative motion 
between the halo neutron and the core nucleus. We assume that the projectile moves along 
a straight line trajectory with a constant velocity. The target nucleus exerts a Coulomb 
potential V^xt(r, t) on the projectile. The potential is treated as a time-dependent external 
field in the projectile rest frame, 

I m„+Mc ' I I ' I 

Here b is the impact parameter which specifies the straight line trajectory, v is the incident 
velocity of the projectile, and Zt and Zc are the charge numbers of the target and core 
nucleus, respectively. The masses of the neutron and the core nucleus are denoted by m„ 
and Mc, respectively. 

The wave function is expanded in partial waves as 

*(r,t)=E^^^>^.n(i-). (3) 

Im 
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When the internal Hamiltonian -^(r) includes the spin-orbit interaction, it is convenient to 
couple the spin of the neutron with the relative orbital angular momentum / to the total 
angular momentum j. Since the generalization to such a case is straightforward, we develop 
the formulation by assuming the wave function of Eq.(3). 

We descretize the time variable in a step At, and represent the wave function of n-th 
time step as *(")(r). The wave function of {n + l)-th time step is calculated by the following 
formula, 

Q-'H{r)At/h^-iV,^t{r,t)At/h^{n)^^-^^ 

The radial part of the wave function of n-th step is denoted as u^l^ (r) . The time development 
is then achieved by two successive procedures. First, the development due to the external 
field Kxt(i", ^) is approximated by 

(-+1)^.^ _,At/;.^(/m|Kxt(r,t)|/W)4':i,(r). (5) 



I'm' 



The evaluation of the matrix element in Eq. (5) is done by expanding the external field into 
multipoles. Next the time development by the internal Hamiltonian is performed separately 
for each angular momentum channel, Im, by using the following approximation 

'"'r"«^i^^.'r*>(.-), (6) 

with 

To obtain i^'i^^^ij') from u']^^\r'), we discretize the radius variable r on mesh points of an 
equal spacing and employ the Crank-Nicolson formula using three-point formula for the 
second order differential operator. 

The time evolution is calculated according to the above prescription from an initial wave 
function, \E'(r,t = — cxd) = ^olr), which is the ground-state wave function of the internal 
Hamiltonian. The breakup component of the wave function is obtained by eliminating all 
the bound state components of the Hamiltonian, 

l<^(r,t)) = (l- E \<l>^)^H) l^(r,t)). (8) 

igbound 

Here the subscript b is kept to stress that the time evolution of the wave function is calculated 
for each impact parameter of the external field. 

The momentum distribution of the relative motion between the neutron and the core 
nucleus after the breakup is obtained by 

^%|^=lim|(k|vI.nr,t))P, (9) 
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where |k) represents the plane wave state. The integration of Eq. (9) over k yields the 
Coulomb breakup probability, Pbu(^) = liiiii^oo(^b^(r, t)|^b^(r, t)), for a given b. Inte- 
grating Pbu(^) over the impact parameter vector yields the total Coulomb breakup cross 
section. 

The breakup cross section can be expressed as a function of the relative motion energy 
E between the neutron and the core nucleus as follows, 



dE 



dhh fdkS(E,-E)^^^^^, (10) 



where E^. = The convergence of the integral in Eq. (10) is very slow with respect to 

the impact parameter. Furthermore the calculation of the relative momentum distribution 
requires a long time step for large impact parameters. To circumvent this difficulty, we 
divide the integration interval of b to two parts, [6mm, and [&s,oo], and in the latter 
interval employ the first order perturbation theory to calculate the energy distribution. The 
value of bs is chosen in such a way that the first order perturbation theory (PT) and the 
time-dependent calculation (TD) give approximately the same energy distribution at 6 = 6^. 
Equation (10) is then recast to 

dE dE dE ' ^ ^ 

The second term of Eq. (11) is expressed in a closed form by using the perturbation theory. 
For this aim we use the breakup probability distribution which is obtained in the first order 
perturbation theory, 

dPSj^\b) WttZt'c'I , dBjEl) 

-^^ = ^Jh^^b^ + (12) 

where Kq and Ki are the modified Bessel functions and ^ — By integrating Eq. (12) 
over the impact parameter in the interval [bs, oo], we obtain 

-^«^^^^'^.6if.R.)^.«.)^. (13) 



dE 9 (;,y)2-^— v..y".v^,*/ 



where ^, = ^. 



hv 

The B{E1) strength function in Eqs. (12) and (13) is defined by 



= ;^|(0,|DlM|0o)|'5(^^-^)+E /^k|(0k|AM|0o)|'^(^k-^) 
d^ Mi M 



M 



--EM^opL p I^. D,m\<I>o) (14) 

TT ^ ti — n -\- IS 



with the dipole operator Dim defined by 
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Dim = — ^ erYiM(r), (15) 
Ac + -L 

where Ac is the mass number of the core nucleus. The wave functions 0, and 0k of Eq. (14) 

are the boTind excited states and continuum states of the Hamiltonian H{y) and Ei and E]^ 
are the corresponding eigenvalues, respectively. To calculate the B[E1) strength function 
we rewrite Eq. (14) in a time- dependent form, 

= |-Re(0o|/^Io/°°cItet(^-^^^)Vi^*D,o|0o) 

= ^Re dtei^'^^^^^^UDl.m)). (16) 

where the ground state is assumed to have / = 0. The wave function il){t) = e~K'f^*Z^io0o 
satisfies the time-dependent Schrodinger equation with the Hamiltonian H{r) and its initial 
wave function ^(0) is equal to -Dio^o- By solving the time-dependent Schrodinger equation 
without an external field, it is possible to obtain the wave function ■0(i) and then calculate 
the B{E1) strength function according to Eq. (16). It is straightforward to generalize Eq. 
(16) to the case where the spin-orbit potential is included. 



III. QUANTUM-MECHANICAL ANALYSIS OF COULOMB BREAKUP 

PROCESS 



To elucidate the breakup mechanism of a neutron-halo nucleus, we investigate the 
breakup process for various choices of the internal Hamiltonian H{r). Before showing cal- 
culated results, we first discuss typical intuitive pictures for the breakup mechanism based 
on the classical mechanics. They include: 

1. Direct breakup mechanism 

2. Resonant breakup mechanism 

3. Free-particle breakup mechanism 

In the direct breakup mechanism, the neutron is assumed to be removed suddenly from the 
core nucleus when the projectile nucleus approaches the point closest to the target nucleus. 
After the breakup occurs, only the core nucleus is accelerated by the target Coulomb field. 
The relative momentum between the neutron and the core nucleus is finite both for the 
longitudinal and transverse directions. When the projectile nucleus has such a resonance 
state that can be excited from the ground state by the target Coulomb field, the resonant 
breakup process may become important and the breakup may proceed dominantly by way 
of the resonance. If the lifetime of the resonance is long enough, the neutron and the core 
nucleus would move together for a long time before the breakup occurs. In such an extreme 
case there is no difference in the relative momentum between the neutron and the core 
nucleus. In the last case of free-particle breakup mechanism we consider the limiting case 
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where the binding energy of the neutron is extremely smalL Then the neutron and the 
core nucleus behave independently during the reaction. The core nucleus moves along the 
Rutherford trajectory, while the halo neutron, receiving no Coulomb force from the target, 
moves on the straight line. In the end the core nucleus receives the momentum only in the 
transverse direction. The momentum difference between the neutron and the core nucleus 
arises only in the transverse direction. Table I summarizes the momentum differences for 
the three cases discussed above. 

It is not clear how well the above classical arguments are quantitatively correct, though 
they are easily accepted intuitively. In what follows we show results of the quantum- 
mechanical calculations for various cases of the neutron-core interactions for which the 
above-mentioned reaction mechanisms are expected to be manifest. 

We choose various parameters for the study of the breakup reaction of ^^Be on a ^ospb 
nucleus performed at the incident energy of 72 MeV/nucleon [§]. Hence the projectile nucleus 
of ^^Be has the incident velocity of f/c = 0.37. By choosing the reaction plane to be x-z 
plane, the target nucleus moves on the straight line RtI^^) = 0, — f t) in the projectile 
rest frame. The impact parameter h is fixed to 6 = 12 fm in the present section. The 
radius variable r is taken up to 800 fm and it is descretized with the mesh size of Ar = 0.4 
fm. The time step At/h = 0.01 MeV~^ is used for calculating the time development of the 
wave function. The Schrodinger equation is solved for the time interval of —10 < t/h < 10 
MeV^^. At the initial and final stages of the calculation, the target nucleus is apart from 
the projectile's center-of-mass coordinate by about 750 fm in the longitudinal direction. 

We assume that the ground state of the ^^Be nucleus is described with a single 
neutron-halo structure in the Is orbit around the inert ^°Be core. The potential between 
the neutron and the core nucleus is taken to be a spherical Woods-Saxon potential. The spin- 
orbit interaction is turned off in the calculation presented in this section. The radius and 
diffuseness parameters of the potential are fixed to R = 2.67 fm and a = 0.6 fm, whereas 
the depth of the potential is treated as a variable parameter to investigate the breakup 
mechanism. The partial waves up to / = 4 are included in the expansion of Eq. (3). The 
contribution of higher partial waves is found to be negligible. The target Coulomb field of 
Eq. (2) is expanded in multipoles around the projectile's center-of-mass, and the dipole and 
quadrupole multipoles are included in the calculation. The dipole field plays a dominant 
role in the present system. The contribution of the quadrupole field is found to be small. 

The average value of the relative momentum between the neutron and the core nucleus 
is calculated by the following formula 

_ (1-r(r.i)|--^V|1-E"(r,t)) 

{<I.BU(r,t)|„BU(,_()) • (17) 

Similarly the average values of the longitudinal and transverse momentum differences, {k\\) 
and are defined by the breakup component \l/^^(r, t). We use the convention that the 
longitudinal and transverse directions indicate the z and x directions, respectively. 
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A. Dependence on the neutron binding energy 



To discuss the validity of the classical arguments summarized in Table 1, wc first show 
the result of calculation obtained by changing the neutron-core potential depth Vq oi I — 
channel. The value of Vq is varied to understand the dependence of the breakup reaction 
mechanism on the binding energy of the halo neutron of the Is orbit. The potential depth 
of other channels, denoted Vp, is set the same as Vq. The Op orbit is always bound below 
the Is orbit, and there is no resonant nor bound excited state which can be excited from 
the ground state by the dipole field. The halo nucleus considered in this subsection is thus 
a very simple system, and excited to the continuum directly by the Coulomb field. 

Figure 1 displays the magnitude of {k\\) and {k±) values as a function of Vq. The lower 
part of Fig. 1 shows the binding energy of the Is orbit. As the binding becomes weaker, the 
longitudinal momentum difference decreases and approaches zero in the vanishing binding 
energy, while the transverse momentum difference increases shghtly. This result is consistent 
with the free-particle breakup picture. However, the transverse relative momentum is always 
smaller than the value, 2kc = 0.092 fm~^, which is expected from the classical argument. 
On the other hand, in the strong binding case, both of the longitudinal and transverse 
momentum differences are finite and close to the value of kc = 0.046 fm~^, which is just 
expected from the direct breakup picture. These calculations including both weak and 
strong binding cases indicate that the breakup mechanism of the neutron-halo nucleus, 
whose binding energy is typically less than 1 MeV, proceeds between the two mechanisms of 
the direct and free-particle breakup. The free-particle picture becomes more suitable with 
the decreasing binding energy of the halo neutron. 

We show in Fig. 2 the B{E1) strength as a function of the neutron-core relative energy. 
The value of e of Eq. (16) is set 0.01 MeV. The B{E1) strength function is closely related to 
the energy spectrum of the breakup cross section (see Eq.(13)). The non-perturbative effect 
gives only a small effect on the energy spectrum of the cross section [14,13b]. The B{E1) 
strength is shown for two cases of the ground-state energy: E = —2.00 MeV (Vq = —65.2 
MeV), and E = -0.503 MeV {Vq = -58.3 MeV). The latter value is chosen to fit the 
empirical neutron separation energy of ^^Be. No resonance exists in p-wave for both cases. 
Therefore the peak in the B{E1) strength which appears at low excitation energy in case of 
E = —0.503 MeV has nothing to do with any resonant character. We see that the energy 
spectrum of the B{E1) strength is very sensitive to the neutron binding energy. Since the 
breakup probability is sensitive to the El strength at low excitation energy, the breakup 
cross section is also sensitive to the neutron binding energy. 



B. Dependence on the level structure of excited states 

Here we investigate how the breakup mechanism depends on the structure of the excited 
states of the halo nucleus. For this purpose we employ the angular momentum (/-) dependent 
neutron-core potential. The potential depth of / = channel is fixed to Vq = —58.3 MeV 
in order to fit the energy of the Is orbit to Eq — —0.503 MeV, the ground-state energy of 
^^Be from the n+^°Be threshold. The value of Vp is now treated as a variable parameter. 
As shown in the lower part of Fig. 3, there is a bound Op state when Vp < —31.3 MeV. In 
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case of —33.2 MeV < Vp < —31.3 MeV, the bound Op state appears between the Is ground 
state and the neutron threshold. This Op state is excited by the dipole field in the breakup 
reaction. In case of —31.3 MeV < Vp < —30 MeV a p-wave resonance with a long lifetime 
appears in the continuum. 

Figure 3 plots the average relative momentum as a function of Vp. As expected, when 
the resonance state exists both of the average relative momenta in the longitudinal and 
transverse directions are very small. It is important to note, however, that the transverse 
momentum remains finite even when the resonance energy is very close to zero and thus its 
lifetime becomes very long. This indicates that the resonant breakup mechanism is too much 
oversimplified. When the potential depth is deep enough to have a bound excited state, 
the average relative momentum shows quite different behaviour between the longitudinal 
and transverse directions. The longitudinal momentum difference is still very small and 
changes continuously from the case where the excited state is the resonance. In contrast 
with this the transverse momentum difference increases discontinuously from the resonance 
case. These results of small longitudinal and large transverse relative momenta might suggest 
that the breakup proceeds through the free-particle breakup mechanism in this case. There 
is, however, no physical reason that the free-particle breakup mechanism is correct. This is 
because two breakup processes occur when there is a bound excited state: One is the direct 
breakup to the continuum. The other is the breakup via the bound excited state. The latter 
process apparently does not fit in with the free-particle breakup mechanism. At present we 
do not have a simple explanation for the discontinuous change in the transverse momentum 
difference and the continuous behaviour in the longitudinal momentum difference which show 
up when the excited Op orbit crosses the threshold. When the Op orbit is bound more deeply 
{Vp < —35 MeV), or is in the nonresonant continuum (Vp > —25 MeV), the momentum 
differences become similar to those of the previous subsection (Vp = Vq = —58.3 MeV, see 
Fig. 1), that is, = 0.026 fm'^ and {k±) = 0.059 fm-\ and are rather insensitive to the 
potential depth. 

The B{E1) strength function is compared in Fig. 4 for two choices of Vp-. One is \^=— 30.0 
MeV where the Op resonance appears at about 0.2 MeV. The other is ^ = —32.0 MeV which 
locates the Op excited state at £'=—0.183 MeV. The latter case corresponds to the previous 
calculation in [Q. This was chosen because the | excited state of ^^Be is known to have 
the strong El transition strength. The B{E1) strength calculated with the /-independent 
potential is already presented in Fig.2. When the resonance exists, the sharp peak appears 
at the resonance energy. Except for the energy region of the sharp peak, the B{E1) strength 
is rather similar between the two cases of the resonance and the bound excited state. The 
B{E1) strength function with the /-independent potential has also similar shape, although 
its magnitude is larger by a constant factor. This difference is understood by noting that, 
for the case of the /-independent potential, there is no resonant nor bound excited state 
in the p-wave so that all the El strength appears in the continuum. The similarity in the 
shape of the B{E1) strength for three cases indicates that the energy dependence of the 
non- resonant part of the B{E1) strength is mainly determined by the wave function of the 
ground state, and is rather insensitive to the level structure of the excited states. 

IV. ANALYSIS OF ^^Be COULOMB BREAKUP REACTION 
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In this section we analyze the Coulomb breakup reaction of ^^Be + "^^^ Pb at the incident 
energy of 72 MeV/nucleon done at RIKEN ^j. To reproduce the known properties of the 
^^Be structure, we choose the potential between the neutron and the ^°Be core nucleus in the 
following way: The depth Vq is determined to reproduce the neutron separation energy of 
the ^^Be (0.503 MeV) by assuming the Is orbit of the halo neutron. For / 7^ channels the 
central and Is potentials are included. The strength of the Is potential is set the standard 
value for the p-shell nucleus, Vis =32.8 MeV-fm^. The strength of the central potential is 
then determined to reproduce the observed ^ excited state located at —0.183 MeV from 
the neutron threshold. The state is assumed to be described simply with the Opi orbit. No 
deformation or clustering effect of the ^°Be core nucleus is taken into account. 

The time evolution of the wave function is calculated by using the same parameter sets 
as the previous section. The Coulomb breakup cross section is obtained by the integration 
over the impact parameter larger than bmm = 12 fm (see Eq. (10)). The reaction of the 
impact parameters smaller than 6min proceeds by the nuclear force as well as the Coulomb 
force, and is assumed to lead to more violent nuclear reaction processes. 

The Coulomb breakup process is dominated by the dipole component of the target 
Coulomb field, and proceeds through the excitation to two angular momentum channels, 
j'^ = ^ and I . The B{E1) transition strength to the bound level of | is calculated to 
be 0.251 eHm^ in the present model, which is considerably larger than the measured value, 
0.115±0.011 e^fm^ |^. The transitions to the other bound levels below the ground state, 
which should not occur in principle, are not excluded in the process of the time evolution of 
the wave function, but the mixing-in of those states is found to be negligible in the present 
calculation. 

We show in Fig. 5 the impact parameter dependence of the breakup cross section 
daBv/db = 2TrbPBv{b). The calculated distribution is compared with the data |^ which 
are extracted from the measured breakup cross sections by using the classical argument 
for the trajectory. Figure 6 compares the energy spectrum of the breakup cross section, 
daBv/dE, with the measurement [Q. The calculation reasonably reproduces the measured 
distribution. Dashed curves in Figs. 5 and 6 show the previous results of [|14| obtained with- 
out use of the Is potential, where the Op orbit was fitted to the energy of the bound excited 
I level. Since the El transition to this excited state, though fairly strong, does not lead 
to the breakup reaction, the breakup probability was underestimated in the previous calcu- 
lation. By the introduction of the Is potential the Op orbit now splits into two levels, Opi 
and Ops, in the present calculation. Since the bound Op 3 orbit is located below the ground 
state and the transition to this state is small, most of the El strength in the pa channel is 
distributed in the continuum. This is the reason why the breakup cross section increased in 
the present model which includes the Is potential. The calculated B{E1) strength function 
shown in Fig. 7 also confirms that the inclusion of the Is potential leads to the increase 
in the strength for the same reason mentioned above. The B{E1) strength is compared 
to the experimental data which are extracted from the breakup cross sections by using the 
method of virtual photon spectra [|l^]. The agreement between theory and experiment is 
rather good, which is expected because the energy spectrum of the breakup cross section 
has already shown reasonable agreement as shown in Fig. 6. We note, however, that the 



10 



present model overestimates the B{E1) strength to the | excited state, and do not know 
what effects the overestimation causes on the breakup mechanism. 

We display in Fig. 8 the longitudinal and the transverse momentum distributions of 
the neutron-core relative motion. The impact parameter is set the smallest value, b = 12 
fm. The longitudinal momentum distribution is obtained by integrating the distribution 
of Eq. (9) with respect to and ky, while the transverse distribution is obtained by the 
integration over ky and kz- The longitudinal momentum distribution slightly shifts to the 
negative direction, whereas the transverse momentum distribution shifts to the positive 
direction. The average value of the momentum of the relative motion is {k\\) = —0.019 
fm~^, and {kj_) = 0.053 fm~^, respectively. The measured difference of the longitudinal 
momentum is about 0.04 fm~^, which is close to the value of kc = 0.046 fm~^ expected 
from the classical picture. Though the calculation reproduces a right order of magnitude, it 
explains only a half of the measured value. 

Finally we consider the breakup mechanism referring to the result of the previous section. 
The mechanism is different in two channels, j'^ = | and | , because they have different 

level structure of the excited states. In the j'^ = | channel there is a bound excited state 
which is very close to the threshold. This situation is similar to the case of = —32.0 MeV 
which we already investigated in Fig. 3 with the Is potential turned off. As expected, the 
momentum difference in this channel, especially in the longitudinal direction, is very small. 
On the other hand, there is no bound excited state in the j'^ = | channel. The situation is 
thus similar to the case of = —45 MeV in Fig. 3. The average longitudinal momentum 
difference increases and turns out to be —0.022 fm~^, about a half of the classical value of kc- 
The breakup reaction of ^^Be is thus considered to occur in these two different mechanisms. 



Our previous treatment |jl^ which did not include the Is potential is equivalent physically 
to including only the former process. It is therefore understandable that we obtained the 
very small longitudinal momentum difference in that case. A realistic choice of the internal 
Hamiltonian, particularly the inclusion of the spin-orbit interaction, is very important for a 
quantitative analysis of the breakup reaction of the ^^Be nucleus. 

V. SUMMARY 



We investigated the Coulomb breakup mechanism of the nuclei with single neutron- 
halo structure, focusing on the mechanism which causes the momentum difference of the 
neutron-core relative motion after the breakup. The breakup process was described in the 
framework of the time-dependent Schrodinger equation by treating the target Coulomb field 
as a time-dependent external potential. 

We first discussed the validity and the limitation of the classical arguments made for 
the post-acceleration effect from the quantum-mechanical viewpoint by investigating the 
dependence of the reaction mechanism on the neutron binding energy and the level structure 
of the halo nucleus. We found that the classical arguments do not have quantitative accuracy 
though it is useful to understand the qualitative aspects of the quantum calculation. A 
behaviour which cannot be understood in the classical picture occurred particularly when 
the bound excited state locates close to the neutron threshold. Small longitudinal and large 
transverse momentum differences were found in this case. 
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Wc analyzed the Coulomb breakup reaction of ^^Be + ^°^Pb by employing a realistic 
potential that describes the relative motion between the neutron and the ^'^Be nucleus. A 
spin-orbit potential was included to describe the splitting of the structure of p states. The 
reaction proceeds through two channels, — | and | . Since the | channel has a bound 
excited state near the threshold but the | channel does not have any bound excited states 
or low-lying resonances, the breakup mechanism is rather different in the two channels. By 
including the neutron-^°Be Is potential we took into account, in the present analysis, the 
difference in the breakup mechanism which is sensitive to the level structure of the nucleus. 
By this improvement the post-acceleration effect was enhanced compared to the previous 
case which neglected the Is potential and moreover the magnitude of the energy spectrum 
of the fragments was in reasonable agreement with the measured values. The momentum 
difference in the longitudinal direction was, however, still underestimated by a factor of two 
compared with experiment. 
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Table I 



The classical estimates of the average longitudinal and transverse momenta, {k\\) and 
{k±), of the relative motion after the projectile nucleus fragments into the neutron and the 
core nucleus by the target Coulomb field. Three reaction mechanisms are classified and 
characterized by the classical momentum kc = \mn/{mn + Mc)]ZTZce'^ /hbv, where Zee 
{Znpe) is the charge of the core (target) nucleus, m„ and Mq arc the masses of the neutron 
and the core nucleus, respectively, v is the incident velocity of the projectile nucleus, and b 
is the impact parameter. The zero momentum for the resonant breakup is the limiting case 
of a long lifetime. 
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Free-particle 
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Figure Captions 



Fig.l. The average values of the longitudinal (solid curve) and the transverse (dashed 
curve) momentum difference between the neutron and the core nucleus as a function of the 
neutron-core potential depth Vq. The impact parameter b is set 12 fm. The energies of the 
bound orbits from the neutron threshold are shown in the lower part. 

Fig. 2. The B{E1) strength as a function of the energy E of the relative motion between 
the neutron and the core nucleus. Two choices of the neutron-core potential depth Vq are 
made to locate the Is ground state at -0.503 MeV {Vq = -58.3 MeV) and -2.00 MeV 
(-65.2 MeV). 

Fig. 3. The average values of the longitudinal (solid curve) and the transverse (dashed 
curve) momentum difference between the neutron and the core nucleus as a function of the 
neutron-core potential depth Vp of I ^ channel. The impact parameter b is set 12 fm. 
The solid curve in the lower part shows the energy of the bound Op orbit from the neutron 
threshold. The energy of the Is orbit is fixed to —0.503 MeV as indicated by dashed line in 
the lower part. 

Fig. 4. The B{E1) strength as a function of the energy E of the relative motion between 
the neutron and the core nucleus. Two choices of the neutron-core potential depth Vp of 
/ 7^ channel are made to locate the Op orbit at -0.183 MeV (V^ = -32.0 MeV) and 0.2 
MeV (—30.0 MeV). The dotted curve indicates the strength reduced to one-tenth of the 
result shown by dashed curve. 

Fig. 5. The Coulomb breakup cross section as a function of the impact parameter b. The 
solid curve is the result of the present model, while the dashed curve is the result with the 
/s-potential turned off |T^. Experimental data are from |p. 

Fig. 6. The Coulomb breakup cross section as a function of the energy E of the relative 
motion between the neutron and the ^°Be nucleus. The solid curve is the result of the 



present model, while the dashed curve is the result with the /s-potential turned off 
Experimental data are from [Q. 

Fig. 7. The B{E1) strength as a function of the energy E of the relative motion between 
the neutron and the core nucleus. Experimental data are from 

Fig. 8. The longitudinal (a) and transverse (b) momentum distributions of the relative 
motion between the neutron and the ^°Be nucleus. The impact parameter b is set 12 fm. 
The solid curve is the result of the present model, while the dashed curve is the result with 
the /s-potential turned off fl^ . 
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